Mon. Not. R. Astron. Soc. 000,[T|-?? (2013) Printed 7 May 2013 (MN WjgL style file v2.2) 

On the development of instability of the black hole-torus systems 
and quasi periodic oscillations 
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ABSTRACT 

We present the numerical study of dynamical instability of a pressure-supported relativistic 
torus, rotating around the black hole with a constant specific angular momentum on a fixed 
space-time background, in case of perturbation by a matter coming from the outer boundary. 
Two dimensional general relativistic hydrodynamical equations are solved at equatorial plane 
using the HRSCS to study the effect of perturbation on the stable systems. We have found that 
the perturbed torus creates an instability which causes the gas falling into the black hole in 
a certain dynamical time. All the models indicate an oscillating torus with certain frequency 
around their instant equilibrium. The dynamic of the accreted torus varies with the size of 
initial stable torus, black hole spin and other variables, such as Mach number, sound speed, 
cusp location of the torus etc., but not their instability. The precessing torus not only effects 
the gravitational radiation, but also generates it. On the other hand, the mass accretion rate 
is slightly proportional to the torus-to-hole mass ratio in the black hole-torus system, but it 
strongly depends on the cusp location of the torus. The cusp located in the equipotential sur- 
faces of the effective potential moves outwards into the torus. The dynamical change of the 
torus increases the mass accretion rate and also triggers a runaway instability. Our numerical 
simulations also show that the oscillating relativistic torus could be used to explain the multi- 
ple peaks observed in the black hole high frequency QPOs, responsible for radiation of X-ray 
observed by different X-ray telescopes. 

Key words: general relativistic hydrodynamics: numerical relativity: black hole: torus: insta- 
bility: quasi-periodic oscillation 



1 INTRODUCTION 

The coalescing of compact binaries and their interactions with 
accretion discs are important issues in astrophysics. The re- 
cent review about the black hole accretion disc theory given by 
lAbramowicz & Fragile! ( 1201 31) explains how the interaction of ac- 
cretion disc with a black hole reveals predictions of emission, sig- 
nature of strong gravity, black hole mass, and spin. The merging 
of BH-BH, BH-NS or NS-NS, or accreting a hot matter rotating 
around compact objects may cause a torus to be formed which 
might be responsible for oscillation and emission of BH-torus sys- 
tems. The hot torus around the black hole is thought to be a mecha- 
nism gather the particles and form a jet. The internal shocks inside 
the jets can create hi gh energy particl es and lead to the emission of 
gamma-ray photonsl lMeszaros 120061) . It is important to understand 
the dynamics, formation, and stability properties of a torus to show 
how the gamma rays are formed. 

The rapid destruction of an accretion disc due to high accre- 
tion rate, which is larger than the Eddington rate, was first sug- 
gested bv lAbramowicz et.alj|1983l) . They pointed out that the cusp 
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located in the equipotential surfaces of the effective potential moves 
outwards into the torus. This physical change in the black hole- 
torus system can increase the mass accretion rate and mass of the 
black hole, and it causes a runaway instability. The instability of 
accreting t orus orbiting around the non-rotating black hole was 
studied bv lMontero et.al.1 yOlOJ), and their numerical simulations 
showed that axisymmetric oscillation of torus was exhibited with- 
out the appearance of the runaway instability. It is also indicated 
that the self gravity of torus did not play an important role for the 
occurring of the instability in a few dyna mical time steps. The run- 
away instability was firstly discovered bv lAbramowicz et.al.Ul983l) 
using a pseudo - Newto nian poten tial and studied in rela tivity by 
iFont & Daignej J2002t) . Recently, iKorobkin etal] I T2012T) . for the 
first time, demonstrated the runaway instability by using the fully 
self-consistent general relativistic hydrodynamical evolution. They 
have confirmed that the runaway instability occurred in a few dy- 
namical times. 

The interacting matter with a black hole or falling matter into 
a black hole is a source of gravitation radiation. Perturbation of the 
torus by matter modifies the fluid around the black hole and shock 
may be formed. The rotating shock can cause a matter to have a 
non-linear behavior and its energy heats the gas. So the electromag- 
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netic flares wou ld be produced and these flay ers continue to occur 
for a long time dSchnittman & Krolik 1120081) . Interacting particles 
with shock waves in the equatorial plane produce a continues emis- 
sion ranging from Ulrta-Viole to X-rays if they freely drive away 
from the black hole-torus system. 

It is known that the central region of the black hole-disc sys- 
tem, which produces the Gamma Ray Burst (GRB), is not sig- 
nificantly bigger than the size of the black hole. So the torus is 
a possible physical mechanism to produce GRB and gravitational 
waves which are used to identify the nature of the inner region that 
may be observed from GRB system. The small perturbations on 
the torus create a hydrodynamical instability which is responsible 
for gravitational waves. Its co ntribution to the grav itational wave 
background was estimated by ICoward et.al.1 (120021) . The gravita- 
tional wave radiation from the cente r of GRB, w hich consists of a 
black hole and torus, was studied bv lLeel I I2004) . It was found that 
the gravitational wave is a natural phenomenon in the wobbling 
system. 

Oscillation properties of the black hole-torus system have 
been studied for numbers of astrophysical systems. Investigating 
the oscillating system can allow us to expose discoseismic classes 
which are dominant in relativistic geometrically thin disc. Some 
numeri cal studies re v ealed t he oscillatory modes o f the perturbed 
torus dZanotti et.all 12003k iRezzolla et.ail 120031 ; iMontero et.ail 
120071) . Oscillatio n properties of the pressure-sup ported torus were 
also examined bv lSchnittman & Rezzollal ( 120061) . They had found 
a strong correlation between intrinsic frequency of the torus and 
extrinsic frequency shown in observed light curve power spectrum. 
The analysis of Quasi-Periodic Oscillations (QPOs) was motivated 
by observed data coming from the galactic center, low-mass X— 
ray binaries (Ivan der Klisll2004|) and high fr equency phenomena, 
testing the strong gravity (IBursa et.al. 1120041) in the innermost re- 
gion of accretion disc around the black hole. 

The Observations with X-ray telescopes hav e uncovered the 
existi n g of QPOs of matte r around the black holes JRemillard et.al. I 
1 19991 : Istrohmaverl 1200 lb . The QPOs can be used to test the 
strong gravity and the properties of the b lack hole since they 
are originated close to the black hole JDonmez et.al.1 1201 ll : 
lAbramowicz & Fragile! 1201 31) . In order to understand these phys- 
ical properties, the perturbed pressure-supported torus is exam- 
ined in a close-binary system. In this paper, our main goal 
is to study the instability of a stable torus on the equato- 
rial plane based on perturbation, which is coming from the 
outer boundary of computational domain. In order to do that 
we solve general relativistic hydrodynamical equations in fixed 
space-time for the non-rotating and rotating black holes at the 
equatorial plane. Thus, we investigate the effects of perturba- 
tion and thermodynamical parameters onto the tori's instabil- 
ity, which causes the appearance of GRBs, in a few dynamical 
times. We have perturbed the to r us in a different way from th e 
one used inlMontero et.all J2007h:ISchnittman & Rezzolhl J2006h; 



IZanotti et.alj d2005» : IZanotti & Rezzollj d2003l) : iFont & Daignei 
1 20021) in which it was performed by applying a displacement (or 
kick) on the torus. 

In this work, we study the effect of perturbation onto the torus- 
black hole system by solving the fully general relativistic hydrody- 
namical equations in a fixed background. The paper is organized as 
follows: in section[2] we briefly explain the formulation, numerical 
setup, boundary and initial condition and introduce the initial setup 
of the relativistic torus. Next, Section [3] describes the numerical 
results. We present the outcomes of interaction of the torus-black 
hole system with matter as a result of perturbation coming from 



the outer boundary. QPOs, due to oscillation of the torus, is also 
presented. Finally, Section [4] concludes our findings. Throughout 
the paper, we use geometrized unit, G — c — 1 and space-time 

signature (-,+,+,+). 



2 EQUATIONS, NUMERICAL SETUPS, BOUNDARY 
AND INITIAL CONDITIONS 

2.1 Equations 

We consider a spherical symmetric, non-self gravitating a torus 
with a perfect fluid equation of state in a hydrostatics equilibrium 
around the non-rotating and rotating black holes. In order to model 
the instability of the torus due to perturbation, we have solved gen- 
eral relativistic hydrodynamical equations using perfect fluid equa- 
tion of state. The general relativistic hydrodynamical equations are 



Vm T^ = 0, Vm^=0, 



(1) 



where T 1 *" = phu^u 1 ' + Pg^ v is a stress-energy tensor for a per- 
fect fluid and J M = pit M is a current density. To preserve and use 
the conservative properties of general relativistic hydrodynamical 
equations, Eq[T|is writte n in a conservati ve form using the 3 + 1 
formalism, and we have I IFont et.alfeOOOl) 



dU dF r m^_ 
dt dr 86 



S, 



(2) 



where U is a conserved variable and, F r and F 1 * are fluxes in radial 
and angular direction at equatorial plane, respectively. These quan- 
tities depend on fluid and th ermodynamical v ariables which ar e 
written explicitly in detail in lDonmezl J2004 : IFont et.alj d2000l) . 
S represents the source term which depends on space-time met- 
ric g**" , Christoffel symbol F" v , lapse function a, determinant of 
three metric jij and stress-energy tensor T*" . The conserved vari- 
ables and fluxes also depend on the velocity of fluid and Lorenz 
factor. The covariant components of three-velocity v % can be de- 
fined in terms of four- velocity u M as v l = w'/( M ') ar, d Lorenz 
factor W = au° = (1 — 7yu , V 3 ) -1 . The explicit representa- 
tion of general relativistic hydrodynamical equat ions and their nu - 
merical solutions at equato rial plane are given in lDonmez I (12004) : 
lYildiran & Donmezl (LJoToh . 



2.2 Numerical Setups and Boundary Conditions 

We have solved the General Relativistic Hydrodynamics (GRH) 
equations at the equatorial plane to model the perturbed torus 
around the black hole. The code used in this paper, whi ch carries 
out th e numerical simulation, was explained in detail bv lDonmezI 
( 120041) . which gives the solutions of the GRH equations using 
Marquina method with MUSCL left and right states. Marquina 
method with MUSCL scheme guarantees higher order accuracy 
and gives better solution at discontinuities mostly seen close to 
the black hole. The pressure of the gas is computed by using the 
standard V law equation of state for a perfect fluid which defines 
how the pressure changes with the rest-mass density and internal 
energy of the gas, P = (T — l)pe. After setting up the stable 
torus with a constant specific angular momentum inside the re- 
gion, ri n < r < r ma x, we have to also define the vacuum, the 
region outside of the n n < r < r max at where all the variables set 
to some negligible values. We have used dynamically unimportant 
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numbers for atmosphere, p a tm = 1CP 12 . The Kerr coordinate is 
used to set up the black hole at the center of computational domain 
using an uniformly spaced grid in r and cf> directions. The r goes 
from 2.8M to 200M for non-rotating, and 1.7 M to 200M for ro- 
tating black holes while (j> varies between and 2-n. Typically, we 
use N r X N$ = 3072 X 256 zones in radial and angular directions 
at the equatorial plane. 

The boundaries must have been correctly resolved to avoid un- 
wanted oscillations. So we set up an outflow boundary condition at 
physical boundaries of computation domain in radial direction. All 
the variables are filled with values using zeroth-order extrapolation. 
While the velocity of matter, v r , must be less than zero at close to 
the black hole, it should be bigger than zero at the outer boundary 
of computation domain. The periodic boundary condition is used 
in angular direction. 



2.3 Initial Torus Dynamics 

The analytic representation of the non-self-gravitatin g relativistic 
torus for a test particle was first examined by lAbramowicz et.al.1 
1 119781) . They found a sharp cusp for marginally stable accreting 
disc which was located at an equatorial plane. Later, the torus with 
a perfect fluid equation o f state was di s cussed in d e tail in Refs. 
IZanotti & Rezzollj ( [2003l) : IZanotti etH] < l2003l . l2005h : lNagar eUri] 
( 20051) . numerically. The matter of torus is rotating at a circular 
orbit with non-geodesic flow between ri„ and r ma x, and the poly- 
tropic equation of state, P — Kp v , is used to build an initial torus 
in hydrodynamical equilibrium with the values of variables given 
in Table [TJ Using the V = 4/3, it mimics a degenerate relativistic 
electron gas. The significant internal pressure in the torus can bal- 
ance the centrifugal and gravitational forces to maintain the system 
in hydrodynamical equilibrium. 

Almost all the astrophysical systems, which have a torus type 
structure, are expected to form as a consequence of the merger of 
two black holes, two neutron stars, their mixed combination or form 
in the core collapse of the massive stars. These phenomena suggest 
that the torus might be in non-equilibrium state after it is formed. In 
order to understand the physical phenomena of a non-stable black 
hole torus system and perturbation of the torus by a matter which 
is coming from outer boundary of computational domain, we have 
considered different initial conditions and perturbations given in 
Table [TJ 



3 NUMERICAL RESULTS 

The numerical results given in this paper are classified depending 
on the torus-to-black hole mass ratio and black hole spin as well 
as the size of the pressure-supported stable torus around the black 
hole. The perturbing torus will uncover its oscillating properties of 
and dynamical exchange during the evolution. 

We adopt the perturbation from the outer boundary at 3.1 < 
4> < 3.18 with parameters, sound speed, and Mach number given 
in TableQJand rest-mass density of perturbation is p v — p c — 0.1p c 
for all models. The perturbation-to-torus mass ratio is M p /Mt ~ 
0.00001. The matter is injected toward the torus-black hole system 
with radial velocities, V r = 0.002, V r = 0.02, V r = 0.2, and 
V r = 0.0001 and with an angular velocity, V* = 0.00001 for all 
V r , accelerated radially and reached the torus around t — 3000A/. 
The ratio of angular velocity of the perturbed matter to Keplerian 
angular velocity is V^/Qk = 1/35 at the outer boundary, r = 
200Af . We use such a small angular velocity for perturbation due to 



the technical difficulties. If the perturbing stream carries significant 
angular momentum, the code crashes. The period of matter at the 
center of torus is around T c = 170Af . The perturbation starts to 
influence the torus after ~ 16 orbital period of the torus which has 
highest density around r ~ 9A/. But it is ~ 4 orbital period for 
models Pi and P2 given in in Table [TJ These orbital periods are 
measured at the location where the rest-mass density is maximum 
in the unperturbed torus. 

3.1 Perturbed Torus Around the Non-Rotating Black Hole 

The perturbations sent into the black hole-torus system trigger the 
matter falling into the black hole due to the instability causing the 
sudden X-ray flare. The observed X — rays which hit the local 
gas of the torus produce light echo. Because the torus around the 
non-rotating and rotating black holes may shrink and come closer 
to the black holes. Hence, it constitutes more hot region to generate 
higher X-ray fluxes. 

To analyze the dynamics of the torus and its instability after 
the perturbation, we plot the density of the torus for the inner region 
at different snapshots for model Pi seen in FigQJ The torus, ini- 
tially stable, is perturbed by matter coming from the outer bound- 
ary. In the color scheme while the red is representing the highest 
density, the blue shows the lowest value of the density of torus. At 
t ~ 3000A/, perturbation hits the pressure-supported stable torus 
and then the spiral structure with a quasi-steady state case is de- 
veloped as time progresses. Once the perturbed torus reaches the 
quasi-steady state, the non-axisymmetric dynamical features are in- 
troduced in the accreted torus. This quasi-steady state structure ro- 
tates around the black hole and produces quasi-periodic oscillation. 
The rotating gas represents wave-like behavior and it travels inward 
or outward. From Fig[TJ we conclude that the torus stays in steady 
state and the inner radius is located around r — 15 M before it in- 
teracts with perturbation. After the perturbation, the inner radius of 
the torus gets closer to the black hole and, at the same time, the cusp 
located in the equipotential surfaces of the effective potential moves 
outwards into the torus. As a result, the inner radius and cusp loca- 
tion of the newly developed quasi-steady-state torus equal to each 



other and oscillate between the points, r c 



4.71M, 



3.48M, 6.2M, 3.63Af and AM which can be seen in Fig.QJand 
the right panel of Fig[2] during the evolution. The specific angular 
momentum corresponding to these cusp locations are t — 3.77 M, 
4.36AT, 3.67A7, 4.23M and AM, respectively. The location of the 
cusp of the perturbed torus moves out from the black hole and 
hence, the accreted torus reaches to a new quasi-equilibrium state. 
In addition to evidence supplied by Fig. [TJ and the right panel of 
FiglU Fig[5]also shows that the angular momentum distribution of 
the torus increases outwards from r ~ 5.25M, which represents 
the location of the cusp at a fixed time, to the larger radii. After 
all, the distribution of density is non-axisymmetric. It is also seen 
in the color plots that the perturbed torus has an oscillatory behav- 
ior, and the oscillation amplitude stays almost constant during the 
compression and expansion of the torus. 

The rotating torus has centrifugal forces which pull out gas 
outwardly around the black and causes less the gas accreted than 
the Bondi rate during the perturbation. The angular velocity of the 
torus is larger than Keplerian one, and it exponentially decays for 
the larger r. The distribution of angular momentum is maintained 
by the disc pressure. The angular momentum of the torus close to 
the horizon may be transferred to the black hole after the torus is 
perturbed. This transformation can be used to explain the spin up 
of the black hole. 
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Table 1. Initial models of the perturbed torus around the black hole at the equatorial plane. From left to right: P shows the model name, a is the spin parameter 
of the black hole, Mt/Mgu is the toras-to-hole mass ratio, T is the adiabatic index, K(geo) is the polytropic constant in geometrized unit, £rj is the constant 
specific angular momentum, r cusp (M) is the cusp location, p c (geo) is the density at the center of the torus, r c (M) represents the location at where density 
is maximum, C' s and M are the speed of sound and the Mach number of perturbation, respectively. Perturbation shows the type of perturbation applied on 
the torus. The total time of the simulations varies depending on models. 



Model -^ 






K(geo) 



Pc(geo) 



io r cusp (M) r c (M) 



Parameters 

Cs(c) 



of Injected Matter 

M Perturbation 



Pi 
Pi 
Ps 

Pi 
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Pi 
Ps 
Pj 



0.0 
0.0 
0.0 
0.0 
0.0 
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Figure 2. Time evolution of rest-mass density along the r for model P3 (left 
panel) and model P\ (right panel). 



In order to reveal the effect of the size of torus and inner ra- 
dius onto the instability and dynamics of the torus, we consider the 
perturbations of the tori for the same mass, but with different sizes. 
The time evolution of the central rest-mass density of torus, each 
of them plotted along r at cf> = 0, are shown for models Pi and P3 
in Fig|2]- For these models, we have applied a perturbation, which 
is p p — p c — O.lpc, from the outer boundary between t = 0M 
and t = 80M to evolute the influence of perturbation to the overall 
dynamics of torus. The perturbation destroys the torus, triggers the 
instability and can cause the location which has the maximum den- 
sity to approach the horizon. Such perturbation induces the mat- 
ter to fall into the black hole or out from computational domain. 
While the maximum rest-mass density of torus slightly oscillates 
for model Pi , seen in the right panel of Fig|2] the rest-mass density 
significantly decreases, seen in the left panel of the same graph, 
during the time evolution. 

The different initial parameters of torus, such as the specific 
angular momentum, the cusp location, the size etc., can cause dif- 



ferent types of dynamics after interaction with a perturbation.The 
perturbed torus given in model Pi forms a new cusp location, 
which is bigger than the cusp location of initial torus. So, the cusp 
located in the equipotential surfaces of the effective potential moves 
outwards into the torus. These dynamical change of the torus in- 
creases the mass accretion rate in a few dynamical time steps, and 
finally, it causes a runaway instability. The computed accretion is 
a result of the decreased specific angular momentum of the torus 
after it is perturbed. The origin of instability might be raised due 
to new values of physical parameters of torus. Thus, the cusp loca- 
tion, the size, and specific angular momentum of torus play a vital 
role in the ons et of instability. These re sults are consistent with the 
work done by JAbramowicz et.al.ll983l) . They have also confirmed 
that the cusp moving outwards into the torus increases the mass 
accretion rate and the runaway instability is triggered due to the 
increased mass of the black hole. In our simulations, we reach the 
same physical conclusion just simply perturbing initially steady- 
state torus instead of changing space-time metric. 

Accretion rate shows time variation in a great amplitude simi- 
lar to one happening in low or hard X-ray time variation. Comput- 
ing the accreting and the mass accretion rate are important indica- 
tors to understand the disc behavior and its instability around black 
holes. The non-axisymmetric perturbation onto the steady- state 
torus, which has a constant specific angular momentum, around 
the black hole can describe the Papaloizou-Pringle instability. In- 
stability of the accreted torus can be subject to perturbation coming 
from the outer boundary in a number of circumstances which may 
cause unstable mass accretion rates ( IPapaloizou & Pringle II 19841 . 
Il985f) . We have confirmed that the accretion rate is governed by 
some physical parameters: the torus-to-hole mass ratio, the loca- 
tion of cusp, initial specific angular momentum of the torus, Mach 
number of the perturbed matter and the angular momentum of the 
black hole. The mass accretion rate computed from Pi, P3, and 
Pi are depicted in the left part of Fig|3] We have noticed that the 
mass accretion rate for the models are the same at early times of 
simulations, but later the models, P3 and P4, drop exponentially as 
a function of time. Interestingly, the mass accretion rate for only 
model Pi does not increase in amplitude, and it oscillates around 
dM/dt = 5xl0" 5 . 
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Figure 1. The rest-mass density of the perturbed torus during the time evolution for Model Pi , given in a linear color scale. The torus is perturbed by a matter 
coming from outer boundary. The matter interacts with the torus around t = 3034M and distorts it. The instability is formed at t = 4540M and stays 
during the evolution. The structure of the torus is highly turbulent and makes quasi-periodic oscillation. The domain is [X m i n , Y m i n ] — > [X max , Y m ax )] = 
[-100M, -100M)] -> [lOOAf, 100A/]. 



The mass accretion rate of the perturbed torus produces lu- 
minosity associated with the central engine of gamma ray burst. 
The expected maximum accretion rate fro m the gamma ray burst 
is as high as dM/dt w 0.01 - lM /s fchen & Beloborodov I 
1 20071) . Fig[3] indicates that the computed mass accretion rate from 
our numerical simulations, 4M. ~ 2 x 10 5 ( — ^ ) ^£f (geo) where 
^jf(geo) is the mass accretion rate in geometrized unit, is in or- 
der of or slightly higher than the expected maximum accretion rate. 
It is also shown in the left panel of Fig[3] that the mass accretion 
rate suddenly increases after perturbation reaches to the torus and 
shortly after the Papaloizou-Pringle instability is developed. All the 
models in Fig[3] exhibit the Papaloizou-Pringle instability which 
grows exponentially with time and present a non-linear growth rate. 
After long time later (it is approximately t — 15000A/ for mod- 



els P5 and Py), it reaches a saturation point and finally reaches a 
new quasi-equilibrium point. Before the saturation point, the torus 
looses mass and then it relaxes to the quasi-stationary accretion 
state. 

It is obvious from the right part of Fig|3]that the accretion rates 
exponentially decay during the time and is always bigger when the 
mass ratio of the torus-to-black hole is larger. Decaying a mass ac- 
cretion rate reaches a constant value dM/dt ~ 0.1Mq/s around 
t — 14000M (t = 82t or b). It is in good agreement with the ex- 
pected accretion rates for accreting disc with the GRBs. As it can 
be seen in Fig|3] the amplitude of mass accretion rate and its behav- 
ior manifest a dependence on the mass ratio of the torus-to-black 
hole and cusp location of the torus. While the mass accretion rates 
for all models decay with time due to the strong gravity which over- 
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Figure 4. The mass accretion rate (left panel) and the rest-mass density of 
torus along r (right panel) are plotted for model P2 . 
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Figure S. The specific angular momentum for models Pi, P§, Pg, and P7 
and for Keplerian are plotted as a function of radial coordinate a long time 
later that the disc is perturbed and instability is developed. 



comes angular momentum of the torus, it oscillates around a certain 
rate for model Pi , shown in Fig[3] The oscillatory behaviors of all 
models are seen obviously during the time evolution. When the ac- 
creted torus gets close to the black hole and cusp moves outwards 
into the torus, the mass accretion rate and instability of the torus 
depend on the cusp location of the torus and the spin of the black 
hole. 

The continuous perturbation of the black hole-torus system 
can lead not only instability of the torus but also oscillates of its 
dynamics around the certain points, seen in Fig|4] The left panel in 
FigS exhibits the time variation of mass accretion rate for model 
P4 computed for r = 8A/ which is true for the perturbed matter 
injecting from the outer boundary continuously. It is seen in this 
plot that the perturbation starts to interact with a torus around t ~ 
3000M , and then nonlinear oscillation is developed. It is noted that 
the continuous perturbation can create different types of dynamical 
structures and accretion rates. 

Accretion into the black hole or out from the computational 
domain is driven by a transport of angular momentum via shock 
waves created on the torus. The higher specific angular momen- 
tum of the torus outwards can allow the torus to become more sta- 
ble. Fig[5]shows the specific angular momentum of the torus along 
the radial coordinate for models P\, P$, Pa, and Py. The Keple- 
rian specific angular momentum is also plotted at the top of these 
models to compare them. The specific angular momentum changes 
drastically depending on the radius. It becomes apparent that the in- 
creasing of the specific angular momentum outwards leads to stabi- 
lize the torus and suppresses the runaway instability. The runaway 
instability visible immediately after the perturbation starts to influ- 
ence of the stable torus. It is also confirmed by Figs[T]and|3] The 
rotating quasi-steady-state torus can only be seen in model Pi. Be- 
cause the substantial amount of the specific angular momentum is 
kept inside the region r < 50A/. Meanwhile, it is seen in Fig[5] 
that specific angular momentum is sub-Keplerian which indicates 
that the gas pressure creates a pressure force and it supplies a radial 
support to the torus. 

The newly developed angular momentum distribution of the 
black hole-torus system, after the non-axisymmetric perturbation, 
might have the local Rayleigh stability condition depending on the 
slopes and its signs. As you can see in Fig(5] the slope of the spe- 
cific angular momentum is bigger than zero (dl/dr > 0), which is 
the limit of the Rayleigh-stable torus, in model Pi for r < 50M, 



but it is dl/dr < and produces the Rayleigh-unstable torus for 
50M < r < 75M. Any increase of the specific angular momen- 
tum in the radial direction stabilizes the torus. 

Fig[6]represents the time evolution of the maximum rest-mass 
density at the center of the torus, each of them normalized to their 
initial values, for models P5, Pa, Py, Pg, and Pg. For these mod- 
els, we have applied the same perturbation with different densities 
(Pp = Pc — 0.1p c ), giving the first kick to the black hole-torus sys- 
tem. The matter begins to influence the torus around i = 3120Af . 
The perturbation triggers the non-stable oscillation of the torus for 
models given in Fig[6] except the model P$. These models expo- 
nentially decay and produce sharp peaks after perturbation reaches 
the torus. Thus, these results imply that subsonic or mildly su- 
personic perturbation produces the pressure-supported oscillating 
torus around the black hole, and amplitude of the oscillation grad- 
ually decreases while the Mach number of perturbation increases. 
Such oscillations can cause the matter to fall into the black hole 
or outward from the computational domain. These physical phe- 
nomena can significantly reduce the rest-mass density of the torus 
during the time evolution seen in Fig(6] Depending on the initial 
configurations, such as, mass ratio of the torus-to-hole and inner 
radius of the last stable orbit of the torus, the reducing rate of rest- 
mass density is defined at the end of simulation. The density of 
the stable initial torus having an initial radius ri„ = 5.49M given 
in model P5 decreases slower than the one with the smaller in- 
ner radius given in models Pa, Py, P$, and P9. This result is also 
confirmed by model Pi . The matter of the torus close to the black 
hole would fall into the black hole due to the strong gravity after 
the torus is triggered by perturbation. Therefore, the mass of the 
black hole and its spin would increase with time significantly. We 
also note that the reduction of the rest-mass density for model Pa 
is slightly slower than the model Py. So the lower the torus-to-hole 
mass ratio looses the faster the matter during the evolution. But it 
needs a further consideration to reveal a global conclusion. 



3.2 Perturbed Torus Around the Rotating Black Hole 

In order to put forth the effect of the rotating black hole onto the 
perturbed torus, we have perturbed the stable torus orbiting around 
the rotating black hole using the initial parameters given in models, 
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Figure 3. Left panel: the mass accretion rates for models Pi, P3, and P4 computed at r = 6Af in geometrized unit. Three different models, the torus with 
different r; n — r c while keeping the same sound velocity or vice verse, are considered to show the effect of perturbation on the torus. Right panel: mass 
accretion rate for models P5, Pg, and P-j. 
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Figure 6. The variation of the maximum density of torus in logarithmic 
scale as a function of time, each of them normalized to its initial value 
p = p c , for the models, P5, Pe, P7, Pg, and Pg. The torus loses the 
matter during the time evolution for all cases while it is almost stable before 
the perturbed matter reaches the torus. 



Pid Piii and P12 in TableQ] We have analyzed the time evolution 
of the rest-mass density of the torus with/without perturbation and 
found that perturbed tori given in models, Pn and P12, become un- 
stable before the perturbation reaches to them but the structure of 
unperturbed torus on a fixed space- time metric background, called a 
Cowling approximation, never changes during the evolution shown 
in FigQ] The top row of Fig(7] clearly shows that the torus in the 
Cowling simulation does not develop a Papaloizou-Pringle insta- 
bility where steady state structure of torus remains almost constant 
during the evolution (i.e. at least t ~ llt or b). Because the neg- 
ligible mass accretion rates through the cusp freeze the growth of 
the Papaloizou- Pringle instability modes in Cowling approximation 
JHawlevll99lh . 

The perturbations on the tori having different inner radii and 
total mass produce diverse dynamics. The small perturbation onto 
the high mass torus, which is more closer to the black hole, in- 
duces the stable torus before the perturbation reaches to it. The spi- 
ral pattern is produced around the rotating black hole, and it causes 
the gas falling into the black hole during the simulations. At the 
same time, it also creates regular oscillation seen in the second and 
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Figure 8. The time evolution of rest-mass density at a fixed point on the 
torus for models, P\\ and P12, for a = 0.9. The inner radii of the models 
are different, but the same perturbation is applied to the both models. 



third rows of Fig(7]and model Pn in Fig[8] It is obvious from the 
model Pi 2 in Fig[8]and, last two rows of Fig|7]that the torus re- 
sponds somewhat differently to the applied perturbation, and the 
torus with a low mass, which is far away from the black hole, has a 
non-oscillatory behavior. Both of these models indicate that at the 
end of the simulation, the rest-mass density of the torus is reduced 
substantially. On the other hand, it is shown that the perturbed torus, 
which has a cusp located inside the black hole, has a non-linear dy- 
namic during the evolution as seen in the last two rows of FigQ] 

3.3 QPOs from the oscillating Torus 

We have performed the simulations using different values of £0, 
r C us P , ri n and r c of the initial torus with a perturbation given in 
TableQ] For these initial perturbed discs, the global oscillations are 
seen at various frequencies. Applying a perturbation to a perfectly 
stable torus gives rise to an epicyclic motion in radial direction due 
to the oscillation mode. The radial epicyclic frequency is the fre- 
quency of the displaced matter oscillating in radial direction. Fig|9] 
consists of a fundamental frequency at ~ 199Hz and number of 
overtones, ~ 399Hz, ~ 600Hz , ~ 800 Hz etc. that they can 
determine the ratios 1:2:3:.... It suggests that the non-linear oscil- 
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Figure 7. Density p of the torus at the equatorial plane for different models, Pin, Pn, and P12 with a black hole rotation parameter, a = 0.9, given a in 
linear color scale. While the top row shows the unperturbed torus around the rotating black hole, the others indicate the perturbation of torus having different 
size and density. The domain is [X min ,Y min ] -> [X max ,Y max )] = [-100M, -100M)] -»• [100M,100M]. 
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lation of the torus due to the perturbation is a consequence of fun- 
damental modes of the torus. The power law distribution given in 
the right panel of Fig.|9]indicates how the peaks appear as a func- 
tion of (j> at fixed r = 8.12A/ which almost represents the point of 
maximum density of the torus. The amplitudes for the correspond- 
ing peaks are almost the same for all cf> at the same frequencies. 
In both graphs, two strong narrow peaks occur at ~ 199Hz and 
~ 399Hz. The QPO behavior of the oscillating torus found in 
our work might be used to explain commensurability frequencies 
observed in t he sources, #1743322, XTE J 1550564, and GRS 
1915 + 105 jMcClintock & Remillardll2004l) . 

The QPO's frequencies from our simulation for the black 
hole M = 2.5 Mp are approximately the s ame as found in 
IZanotti & Rezzollj d2003h : IZanotti et.all J2003I) . They computed 
the power spectrum of the L-2 norm error of the rest-mass den- 
sity for different r\ after the long time evolution. The fundamen- 
tal frequencies, which are seen around 200Hz, and the series of 
overtones for both cases are observed. Our study indicates that 
the perturbed torus, which is caused either by the introduction of 
a s uitable parametrize d pertu r bation to the radial velocity done 
by IZanotti & Rezzollal J2003h : IZanotti et.all J2003h or by a mat- 
ter coming from outer boundary toward the torus (our simulations), 
produces regular oscillatory behavior. The modes observed from 
both perturbations are called the pressure (p) mode of oscillation 
of the torus. The small size of accretion torus orbiting around the 
black hole can be us ed to explain p mode in high freq uency quasi- 
periodic oscillations JAbramowicz & Kluzniakl2004T) . The p mode 
represents the oscillation of matter in the horizontal direction at 
the equatorial plane and is trapped inside the torus. Due to the 
two-dimensional structure of the system, the magnetic fi eld did not 
create a strong effect on the wave properties of p mode JFu & Lai I 
l2009h . 

It is known from previous discussion that the torus loses mat- 
ter during time evolution, and it has chaotic behavior seen in models 
Ps and Pi. Due to this irregular non-linear oscillation, only a fun- 
damental mode appears, but their overtones are absent in the power 
spectrum when the torus is initially more close to the black hole. 
We suggest that the variation of the torus's matter results from the 
non-spherical symmetric perturbation of the inner accreted torus. 

We have also analyzed the time evolution of mass accretion 
rate for the model, Pn reported in Table [T] to compute Fourier 
spectra of the perturbed torus around the rotating black hole. Af- 
ter the perturbation starts to influence the torus dynamics, mass ac- 
cretion rate oscillates around its quasi-equilibrium point. Due to a 
distinctive quasi-periodic oscillation of the unstable torus, the pe- 
riodic character of the spiral structure becomes more evident after 
perturbation reaches to the torus a long time later. After the torus 
relaxes to regular oscillation, it loses the mass, and the oscillation 
amplitude decreases during the evolution. The computing QPO fre- 
quencies from the mass accretion rate might not give definite infor- 
mation, but we may approximately predict it. The power spectra 
of mass accretion rate shows a fundamental mode with a strong 
amplitude at 500Hz and series of overtones located at 1001Hz, 
1491Hz, and 1998Hz for a = 0.9 and M = 2.5 M©. It is worth 
emphasizing that a fundamental frequency from the torus around 
the rotating black hole is much larger than the non-rotating case 
seen in Table [2] This is an acceptable result because the torus is 
more closer, and gravity is more stronger in case of the rotating 
black hole. Table[2]also suggests that a linear scaling is present for 
the black hole spin. 

The observed high frequency QP Os around the black holes 
and neutron stars were described by iKluzniak & Abramowicz I 



d200lh: lAbramowicz&Kluzniakl teOOll): iBlaes et.alj J2007l) : 
iMukhopadhvav I 120091) ; iMukhopadhvav et.alj J2012I) using a sin- 
gle model which has addressed the variation of QPO frequencies. 
Based on the proposed model, they had predicted the spin of the 
black hole using the observed QPO pairs which seem to appear at 
a 3 : 2 ratio. The lower and higher frequencies of pairs for any sys- 
tem are vi = u± — v s /2 and vh = v r + v s , respectively. Here, u±, 
v T a nd v s are theoret ical values of radial, vertical epicyclic frequen- 
cies JDonmez 1120070 and spin frequency of the black hole, respec- 
tively. Using the above model we find that the computed commen- 
surable frequency, 600 : 399, is exposed at ~ 8.1M. The estimated 
value of location of resonance from numerical computation, which 
is ~ 9M, is in the suitable range of theoretical calculation. 



4 CONCLUSION 

We have performed the numerical simulation of the perturbed ac- 
creted torus around the black hole and uncovered the oscillation 
properties to compute QPOs. Instead of giving any kick to a sta- 
ble torus studied by different authors in literature, the matter is sent 
from the outer boundary toward the torus-black hole system to per- 
turb it. We have initially considered the different sizes of the stable 
tori and torus-to-hole mass ratios to expose oscillations of the tori, 
and computed the fundamental modes and their overtones which 
results from the excitation of frequencies due to perturbation. 

We have confirmed that the torus around the black hole might 
have quasi-periodic oscillation only if we choose a suitable initial 
data for a stable accretion torus. There are a number of different 
physical parameters which may affect the oscillation properties of 
the torus. It is seen from our numerical simulations that the size 
of torus affects the frequencies of oscillating torus and their over- 
tones. The size of torus and its rest-mass density also influence the 
power of oscillations and this oscillation creates radiation o f hotter 
photon s. The similar conclusion has been also confirmed by [Bursa] 
( 120051) . The oscillation and mass losing rate of the perturbed torus 
strongly depend on the Mach number of perturbing matter. Upon 
the given appropriate parameters of perturbation, the subsonic or 
mildly supersonic flow produces the pressure-supported oscillating 
torus around the black hole and the amplitude of oscillation grad- 
ually decreases while the Mach number of perturbation increases. 
For the most of the models given in Table [T] we have found that 
the rest-mass density of the torus substantially decreases due to the 
cusp location which is moving outwards into the torus during the 
evolution. Thus, the inner radius, the specific angular momentum 
of torus and the cusp location play a critical role in the onset of the 
runaway instability. The instability can cause the tours to fall en- 
tirely into the black hole. Hence, the mass of the black hole would 
grow. 

Our studies also indicate that QPOs are common phenomena 
on the disc around the black holes. If the accretion disc or torus 
have quasi-periodic behavior, it emits continuous radiation during 
the computation. The amplitude of oscillation is excited by non- 
linear physical phenomena. It has been explained in terms of the 
excitation of pressure gradients on the torus and it is called p mode. 
It is shown that the power law distribution of oscillating stable torus 
includes a fundamental frequency at ~ 199Hz and their number of 
overtones, ~ 399-ffz, ~ QQQHz, and ~ 800Hz that they can de- 
termine the ratios 1:2:3:4. These frequencies are observed at 
any radial and angular directions of the torus. The computed QPO 
freque n cies are almost the same as those of IZanotti & Rezzollal 
( 120030 : IZanotti et.alj ( 120031) even though they apply a kick to the 
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Figure 9. Power spectrum of the rest-mass density at a single point (any <j> or r) for M = 2.5 Mq for model Pi shown in Table[T] A fundamental frequency 
appears at far left of the both graphs. The others are the consequence of overtones. 

Table 2. The frequencies of genuine mode and its series of overtones, due to the nonlinear coupling, inside the torus from the evolution of mass accretion rate. 
It is computed for different black hole spins while the mass of the black hole is assumed to be M = 2.5 Mq . 
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radial velocity of the torus in order to have an oscillation. The mode 
called as p is also the same in both simulations. We have also con- 
firmed that the fundamental QPO frequencies and their overtones 
for rotating black hole are much higher than the non-rotating case 
due to strong gravity and location of the inner radius of the ac- 
creted torus. The strong gravity plays a dominant role in the high 
frequency modulat ion of observed X-ray flux in the binary system 
JBursa et.al. 1120041) . It is consistent with the computation from our 
simulations seen in Table[2] 

We have also reported that the initial perturbation induces the 
torus-black hole system to start the instability around the rotating 
black hole before the transferred matter reaches to the black hole. 
The perturbation can also lead the matter of the torus to fall into 
the black hole. It reduces significantly the total rest-mass density 
of the torus with time, and the oscillatory behavior of the rest-mass 
density changes depending on the inner radius of the torus. 

Finally, we have performed the numerical simulation using 
2D code on the equatorial plane. Clearly, to investigate the effects 
of forces in the vertical direction on the torus instability that may 
play an important role on the dynamics of the whole system, 3D 
numerical simulations are required. The mode coupling in the ver- 
tical direction is likely to affect the instability of a perturbed torus- 
black hole system. Therefore, we plan to model the same perturbed 
system by using 3D code in the future. 
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